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S^ ! Abstract 

The shell correction is proposed in the improved isospin dependent quantum molecular 
dynamics (ImlQMD) model, which plays an important role in heavy-ion fusion reactions near 
Coulomb barrier. By using the ImlQMD model, the static and dynamical fusion barriers, 
t^^ ' dynamical barrier distribution in the fusion reactions are analyzed systematically. The fusion 

and capture excitation functions for a series of reaction systems are calculated and compared 

f~^ ' with experimental data. It is found that the fusion cross sections for neutron-rich systems 

O ' 

increase obviously, and the strong shell effects of two colliding nuclei result in a decrease of 



the fusion cross sections at the sub-barrier energies. The lowering of the dynamical fusion 
j^ ■ barriers favors the enhancement of the sub-barrier fusion cross sections, which is related to 

the nucleon transfer and the neck formation in the fusion reactions. 

PACS: 25.60.Pj, 25.70.Jj, 24.10.-i 

1 Introduction 

Heavy-ion fusion dynamics at energies near and below the Coulomb barrier has been an important 
subject in nuclear physics for more than 20 years [IJ, which is involved in not only exploring several 
fundamental problems such as quantum tunneling etc, also investigating nuclear physics itself, such 
as nuclear structure, synthesis of superheavy nuclei etc. The experimental fusion excitation func- 
tions can be well reproduced by the various coupled channel methods, which include the couplings 



of the relative motion to the nuclear shape deformations, vibrations, and nucleon-transfer. How- 
ever, the coupled channel models still have some difficulties in describing the fusion cross sections 
for very heavy symmetric systems, in which the neck formation will play an important role. The 
microscopic dynamical description of the fusion reactions is very necessary. Microscopic transport 
theories such as isospin dependent quantum molecular dynamics (IQMD) based on QMD model 
[2] or isospin dependent Boltzmann Uehling Uhlenbeck (IBUU) model based on BUU theory [3] 
are suitable for describing the dynamical process of the fusion reactions, which have been used 
successfully to investigate heavy- ion collisions at intermediate energies |H El E]- However, these 
models have some difficulties for studying the fusion reactions near Coulomb barrier, because of 
some unphysical nucleon emissions in the simulation process of projectile and target approaching. 
Wang et al. had made an important improvement to construct stably initial nucleus in the QMD 
model [7], in which the interaction potential was derived from the usual Skyrme interaction and 
the parameters were adjusted by fitting experimental root-mean-square radii, average binding en- 
ergies etc. Recently, fermionic molecular dynamics (FMD) model has also been used to investigate 
the fusion reactions of oxygen isotopes [S] . 

Based on IQMD model, the interaction potential, nucleon's fermionic nature and two-body 
collision have been improved systematically in ImlQMD model [TJ [9]. In this paper, the shell cor- 
rection is further considered, and its influence on the fusion cross sections is analyzed and compared 
with available experimental data. We present a dynamical analysis on the static and dynamical 
fusion barriers, the dynamical barrier distribution, and the fusion (light and intermediate systems) 
and capture (heavy systems) excitation functions and compare them with experimental data. 

In Sec. 2 we give a description on the ImlQMD model. Calculated results of the fusion 
dynamics and the fusion cross sections are given in Sec. 3. In Sec. 4 conclusions are discussed. 

2 Model description 

The same as the QMD or IQMD model, the wave function for each nucleon in ImlQMD is repre- 
sented by a Gaussian wave packet 
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Here the rj(t), pj(t) are the centers of the ith nucleon in the coordinate and momentum space, 
respectively. The L is the square of the Gaussian wave packet width, which depends on the size of 
nucleus. The total N-body wave function is assumed as the direct product of the coherent states. 



in which the anti-symmetrization is neglected. After performing Wigner transformation for Eq. 
(1), we get the Wigner density as 
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The density distributions in the coordinate and the momentum space are given by 



P(r,^) = / f{r,p,t)dp = Y^ 
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respectively, where the sum runs over all nucleons in the systems. Here we have considered the 
uncertainty relation. 

The time evolutions of the nucleons in the system under the self-consistently generated mean- 
field are governed by Hamiltonian equations of motion, which are derived from the time dependent 
variational principle and read [21 HO] 



dH 
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The total Hamiltonian H consists of the kinetic energy, the effective interaction potential and the 
shell correction part as follows: 

H = T + U,nt + Ush. (7) 

For the kinetic energy we can get it from the total wave function 
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Here the first term is the classical kinetic energy. The second term arises from the Gaussian width 
in momentum space, which is usually omitted in QMD or IQMD calculation. It has the value 7.8 
MeV if we take mi=938 MeV/c^ and L=2 fm^. 

The effective interaction potential is composed of the Coulomb interaction and the local inter- 
action 

Uint = Ucoul + Uioc- (9) 

The Coulomb interaction potential is written as 
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where the tzi is the zth component of the isospin degree of freedom for the ith nucleon, which is 
equal to -1 and 1 for proton and neutron, respectively. The Vij = |ri — rj| is the relative distance 
of two nucleon. The second term on the right side is the exchange term with R = (r^ + rj)/2 only 
for protons, which is important for light nucleus. 

In the ImlQMD model, the local interaction potential is derived directly from the Skyrme 
energy- density functional [HI [12], which is expressed as 

Uioc= fVioc{p{r))dr. (11) 



The local potential energy- density functional reads 



^/oc(p) = f ^+T^^+f^(Vp)^+f^[v(p„-p,)]^+(V+ V+^+c^//^)5^+^.//Vpf , 

2 po 1 + 7 Po 2po 2po 

(12) 

where the p„, pp and p = Pn + Pp are the neutron, proton and total densities, respectively, and 

the 6 = {pn — Pp)/{pn + Pp) is the isospin asymmetry. The first two terms are the bulk energy 

term, which are used in the usual QMD or IQMD model. The surface term, surface symmetry 

term and bulk symmetry term are from the third to the fifth term in turn. The last term is the 

effective mass term, which is related to the effective mass of nucleon together with the third part 

of the fifth term. Comparing with the standard Skyrme interaction pTl [T2] , the spin-orbit term 

is neglect in the ImlQMD model, which gives the shell correction in the binding energies. We 

will use a phenomenological expression to embody the shell effect in the ImlQMD model. The 

coefficients in Eq. (12) are related to the Skyrme parameters as 

a 3 (3 ts 

2 = 8^0, YT^ = Y^p2, (13) 

^ = ^i^tl-^t2-^X2t2)po, (14) 

2 d4 

^ = -^[3ti(2xi + 1) + t2(2x2 + l)]po, (15) 

9r = -(3ti + 5t2 + 4x2t2) (-vr^j pf . (16) 

The parameters in the symmetry energy term are given by 

^ = -^(2xo + l)to, 5 = -^(2x3 + l)t3, (17) 

1/3 \^/'^ 
^ = "12 (2''') '^^^^''^ + ^^ " ^'^^"'^ + ^^^- ^^^^ 

The to,ti,t2,t3 and Xo,a;i,X2,X3 are the parameters of Skyrme force. In this work, we only 
consider the first term in the bulk symmetry energy density in Eq. (12) and take the parameter 



Table 1: ImlQMD parameters and properties of symmetric nuclear matter for Skyrme effective 
interactions 



Force 


SkM* 


SIII 


SkP 


RATP 


SLy4 


SLy6 


a (MeV) 


-317.4 


-122.7 


-356.2 


-259.2 


-298.7 


-295.7 


/3 (MeV) 


249.0 


55.2 


303.0 


176.9 


220.0 


216.7 


7 


7/6 


2 


7/6 


1.2 


7/6 


7/6 


^snr(MeV fm^) 


21.8 


18.3 


19.5 


25.6 


24.6 


22.9 


^r(MeV fm^) 


-5.5 


-4.9 


-11.3 


0.0 


-5.0 


-2.7 


9r (MeV) 


5.9 


6.4 


0.0 


11.0 


9.7 


9.9 


as (MeV) 


30.1 


28.2 


30.9 


29.3 


32.0 


32.0 


Poo (fm-3) 


0.16 


0.145 


0.162 


0.16 


0.16 


0.16 


m*^/m 


0.79 


0.76 


1.00 


0.67 


0.70 


0.69 


K^ (MeV) 


216 


355 


200 


239 


230 


230 



A = as/(2po), which corresponds to the form of the linear dependence of the symmetry energy. In 
Table 1 we list the ImlQMD parameters related to several typical Skyrme forces. In the calculation 
we take the Skyrme parameter SLy6, which can give the good properties from finite nucleus to 
neutron star [12j. 

The shell correction energy in the ImlQMD is written as a Woods-Saxon form: 

1 + xJ^ Esh 
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"^ N ^l + exp[(r,-i?o)/a]' 
where the shell correction energy of ground-state nucleus Egh is given by the method used in 



Droplet model [14]. The r, is the radial coordinate in the center of mass system for the ith 
nucleon. The N,Ro,a are the total number of nucleons, nucleus radius and dispersion width, 
respectively, which has the relation Rq = 1.16A^^/'^. The coefficient x is a correction factor coming 
from the Woods-Saxon form, which is determined by 

1 Tf dr a , , 

X = / ^ 1 62 (20) 

RqJt, 1 + exp[(r - i?o)/a] ' Rq' 

We take a = 0.65 fm in the calculation. The Tj and rj are the radial positions at 0.9po and O.lpo 
as shown in Fig. 1 (a) for nucleus ^''^Pb with the Fermi distribution. In Fig. 1 (b) we give a 
comparison of the radial density distributions of the Fermi form, rectangular, SHF with SKM* 
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Figure 1: Comparisons of the radial density distributions of the Fermi form, rectangular, Skyrme 
Hartree-Fock (SHF) with SkM* parameters and the time evolutions in ImlQMD model. 

force and the time evolutions in the ImlQMD model. Combining Eq.(6) and Eq.(19), the shell 
correction slightly constrains the nucleon motion around the surface of a nucleus, which affects 
the fusion dynamics at near Coulomb barrier energies. 

In the dynamical process of two nuclei approaching, we deal with the shell correction of the 
colliding system with a switch function method [Qj [15] which is used in the classical molecular 
dynamics simulation [16]. The shell correction energy of the system in the dynamical evolution is 
written as 
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where the labels z, j denote the sum over the colliding nuclei and the nucleons of the considering 
nucleus, respectively. The factors El^ and Rq are given by 



El^{R)=E,S{R) + E,{l-S{R)) 



(22) 



and 



Rq = RiS{R) + i?e(l - S{R)) 



(23) 



respectively. Here the Ei^c and Ri^c are the shell correction energies and the radii for projectile 
[i = 1), target {i = 2) and compound nucleus (c), respectively. The switch function is expressed 
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Figure 2: The time evolutions of the shell correction energies for the double magic nuclei '^^Ca 
and ^°^Pb (a), the distance of centers between projectile and target (b) and the shell correction 
energies of projectile, target and system in the reaction ^^Ca+^°'^Pb— ^ ^^^No (c) at the incident 
center-of-mass (cm.) energy 200 MeV. 



as 



S = Co + C^x + C2X' + Csx' + C,x' + C,x\ (24) 

where x = ( r .^r ) ^^^ R is the distance of the centers between projectile and target. The 
quantities Rd and Rcf are the distances at the initial time and at the compound nucleus formation. 
In the calculation, we set Rd = 10 + -Ri + -R2 f^i and Rcf = Re — R2 fm. The coefficients Cq, 
Ci, C2, C3, C4 and C5 are taken to be 0, 0, 0, 10, -15 and 6, respectively. Fig. 2 shows the time 
evolutions of the shell correction energy for the double magic nuclei ^^Ca and ^osp]^ using Eq. 
(19) and in the reaction ^^Ca+^°^Pb^ ^^^No using Eq. (21). For the single nucleus evolution 
in ImlQMD model, the calculated shell correction energies approach the empirical formula [2] 
(-5.48 MeV for ^^Ca and -8.93 MeV for ^°^Pb). In the approaching process of two colliding nuclei, 
the shell correction energy of the system evolves from the sum of projectile and target to the one 
of the compound nucleus. 

The phase space constraint method is also introduced in the ImlQMD model, which is proposed 
by Papa et al. in order to improve the nucleon's fermionic nature [T7]. The one-body occupation 



probability is given by 

7i = J2 ^s^s^Ktj / , /j(r, p, t)t/rc?p, (25) 

where the quantities Sj and Tj represent the quantum numbers of spin and isospin for the ith 
nucleon, respectively. The integration is performed within the phase space volume h^ at the 
center position ri, pi. In the dynamical evolution, we timely examine the occupation probability 
with the condition /^ < 1. Otherwise, a series of two-body scattering are performed to reduce the 
phase space occupation probability. 

The main difficulties of the usual QMD and IQMD model in the description of heavy-ion 
fusion reactions near Coulomb barrier are to construct a stable nucleus over the reaction time of 
the compound nucleus formation (typically 10^^^ ~ 10^^*^ s). For that Maruyama et al. proposed 
a cooling method in the study of the fusion reaction ^^O+^^O [18j. The stability of a single nucleus 
is obviously improved in the ImlQMD model. We give an examination of the root-mean-square 
radii and the average binding energies for the magic nuclei as shown in Fig. 3. The stabilities 
maintain 800 fm/c for the selected nuclei. It is necessary for investigating the fusion dynamics of 
light and intermediate nuclei or the capture process of heavy colliding systems. For the complete 
fusion of heavy systems, especially the synthesis of superheavy nuclei, or the damped reactions 
of two very heavy colliding nuclei |T9j, the stability of the single nucleus needs to be keep longer 
time. 

3 Results and discussions 

In this section the ImlQMD model is applied to analyze the fusion dynamics at energies near 
Coulomb barrier. It is well known that the enhancement of the sub-barrier fusion cross sections is 
related to the various couplings of the relative motion to the nuclear shape vibrations, deformations 
and nuclear transfer degrees of freedom. The fusion barrier at various energies and the dynamical 
barrier distribution are investigated systematically using the ImlQMD model. The influence of the 
shell effect on the fusion cross sections is analyzed. The fusion and capture excitation functions 
for a series of reaction systems are calculated and compared them with experimental data. 
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Figure 3: The time evolutions of the root-mean-square radii and the average binding energies for 
the magic nuclei ^^O, ^^Ca, ^ozr, i^sgn and ^ospb. 



3.1 Fusion barrier and dynamical barrier distribution 

The interaction potential V{R) of two colliding nuclei as a function of the distance R between 
their centers is defined as [20] 

V{R) = Ept{R)~Ep-Et. (26) 

Here the Ept, Ep and Et are the total energies of the whole system, projectile and target, re- 
spectively. The total energy is the sum of the kinetic energy, the effective potential energy and 
the shell correction energy. In the calculation, the Thomas-Fermi approximation is adopted for 
evaluating the kinetic energy P]. The static and dynamical interaction potentials are calculated 
according to Eq. (26) as shown in Fig. 4 for head on collisions of the reaction systems ^^Ca+^^^U 
and ^^S+^'^'^^Zr. The static interaction potential means that the density distribution of projectile 
and target is always assumed to be the same as that at initial time, which is a diabatic process 
and depends on the collision orientations and the mass asymmetry of the reaction systems. The 
total density of the system is taken as the sum of the individual nucleus in the overlapping region. 
For a comparison, the proximity results [2T] and the adiabatic barrier as mentioned in Ref. [22] 
are also shown in the figure (left panel) and the corresponding barrier heights are indicated for 
the various cases. However, for a realistic heavy-ion collision, the density distribution of the whole 
system will evolve with the reaction time, which significantly depends on the incident energy and 
impact parameter for a given system. In the calculation of the dynamical potentials, we only 
pay attention to the fusion events as shown later, which give the fusion dynamical barrier. At 
the same time, stochastic rotation is performed for different simulation events. One can see that 
the heights of the dynamical barriers are reduced gradually with decreasing the incident energy 
and increasing the neutron number of the target. The lowering of the dynamical fusion barrier 
is in favor of the enhancement of the sub-barrier fusion cross sections, which can give a little of 
information that the cold fusion reactions are also suitable to produce superheavy nuclei. We can 
understand the microscopic process of the lowering from the neck formation which derives from 
the nucleon transfer and the dynamical deformation in the fusion reactions |23j . 

The dynamical fusion barrier is calculated by averaging the fusion events at a given energy and 
an impact parameter. To explore more information on the fusion dynamics, we also investigate 
the dynamical barrier distribution as shown in Fig. 5 for head on collisions of the reaction 
36g_j_902j^ g^^ incident energies 80 MeV and 85 MeV (static barrier Vb = 84.37 MeV), respectively. 
The distribution moves towards the lower barrier region with decreasing the incident energy. A 
number of fusion events are located at the sub-barrier region, which is favorable to sub-barrier 
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Figure 4: Comparison of the static fusion barriers (including pole-to-pole and waist-to-waist col- 
lisions), adiabatic barrier, proximity results and dynamical barriers at different incident cm. 
energies for the reaction ^^Ca-|-^^^U (left panel), as well as system dependence of the static and 
dynamical fusion barriers (right panel). 

fusion reactions. Rowley et al. proposed a method of extracting the fusion barrier distribution from 
the second energy derivative of the fusion excitation function [24j. A portion of the distribution 
is located at the sub-barrier region from the analysis of experimental fusion excitation functions, 
which may help us understanding the sub-barrier fusion and the kind of various couplings. 



3.2 Fusion and capture excitation functions 

After constructing the stable events of projectile and target, the simulation of the fusion reaction 
can be performed. Stochastic rotation around their centers of mass for each nucleus by a Euler 
angle is made for every event at a given incident energy and an impact parameter. The simulation 
events are set to be 200 for each incident energy E and impact parameter b, and be 300 at the 
sub-barrier energies. The fusion cross section is calculated by the formula 



afus{E) = 271 / bpf^,{E, h)dh = 27iY^ bpfus{E, b)Ab, 

^0 b=Ab 



(27) 



where pfus{E,b) stands for the fusion probability and is given by the ratio of the fusion events 
Nfus to the total events Ntot- The reliability of the calculated pfus{E, b) is estimated by the value 
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Figure 5: Dynamical barrier distribution for the reaction ^^S+^'^Zr at incident energies 80 MeV 
and 85 MeV, respectively. Arrows show the static fusion barriers. 

In the calculation the step of the impact parameter is set to be Ah = 0.5 fm. We also use the same 
procedure in the definition of the fusion event in Ref. [7j that the coalesced one-body density can 
survive through one or more rotations of the composite system or through several oscillations of 
its radius. The capture cross section of heavy colliding system is also calculated using Eq. (27), 
in which the capture event is defined as that the distance between the centers of colliding nuclei 
is smaller than the minimum position of the pocket of the static interaction potential. In Fig. 6 
we show the dependence of the fusion probability as a function of the impact parameter on the 
incident energies in the reaction ^^S+^°Zr and on the system size at the energy £'c.m.=80 MeV. 
One can see that the higher incident energy and the neutron-rich system have the larger fusion 
probability, and the fusion probability is reduced with increasing the impact parameter or the 
relative angular momentum at each incident energy. There is no signal on the appearance of the 
so-called "fusion window" as predicted by the time dependent Hartree-Fock method j25j . 

Recent experimental data showed that there was no enhancement of the fusion cross sections 
in the reactions ^^Ca-|-^°Zr or ^^Zr comparing those in the reactions ^''Ca-|-^°Zr or ^^Zr although 
using more neutron-rich projectile ^^Ca [26]. Possible reason is explained from the fact that the 
double magic nucleus ^^Ca has more rigid (stronger shell effect) structure than ^''Ca. We analyzed 
the influence of the shell effect on the fusion cross sections in the reaction ^^Ca-|-^°Zr as shown 
in Fig. 7. With the same initial conditions such as the stable event and the distance of the 
centers between two nuclei etc, the calculated fusion cross sections are obviously reduced after 
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Figure 6: The fusion probability as a function of the impact parameter and its dependence on the 
incident energies and on the system size. 

considering the shell correction, especially in the sub-barrier domain. The main reason is that the 
shell correction potential in Eq. (7) constrains the motion of the surface nucleons of the colliding 
nuclei. The experimental data can be reproduced rather well by considering the shell correction. 
Fig. 8 shows a comparison of the calculated fusion excitation functions by the ImlQMD model 
after considering the shell correction, one- dimensional Hill- Wheeler formula [27] and the exper- 
imental results for the reactions 46'jj_|_46rpj jjgj g^j^^ ^'^Ca+^^^Sn [29j. The Hill- Wheeler formula 
underestimates the fusion cross sections at the sub-barrier energies for the two reaction systems. 
There is no other adjustable parameters in the ImlQMD model, which is a purely dynamical 
process. The agreement of the calculated fusion excitation functions with the experimental data 
is remarkably well within statistical error bars. For systematically examining the reliability of the 
model and exploring the fusion dynamics, we calculated the fusion cross sections of a series of 
the reaction systems and compared them with available experimental data [261 EOl EH |32l |33] 
as shown in Fig. 9. One can see that the neutron-rich combinations have the larger fusion cross 
section, especially at the sub-barrier regions. As we discussed in the section 3.1, the neutron-rich 
systems give the lower dynamical fusion barriers at the same incident cm. energy, which are fa- 
vorable to increase the fusion cross sections, especially at the sub-barrier energies. The reactions 
i8Q_j_58]\j-j ^^^^ ^^0-|-^''Ni lead to the same compound nucleus formation. However, the system 
i8Q_j_58]\j-j ]-^g^g ^]-^g larger fusion cross sections at the sub-barrier energies. Zagrebaev explained the 
enhancement owing to the positive Q value of the neutron transfer by using a simplified model 
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Figure 7: The calculated fusion cross sections with and without considering shell effect in the 
reaction ^'^Ca+^^Zr and compared with experimental results [26] . 



j. Therefore, the colliding system has a little surplus energy to pass over the interaction bar- 
rier. In the ImlQMD model, since the magic nucleus ^^O has the stronger shell effect than ^^O, 
the transfer of the surface nucleon is slightly constrained in the course of projectile and target 
approaching, which also results in the decrease of the sub-barrier fusion cross sections. 

The synthesis of heavy or superheavy nuclei is mainly reached by the complete fusion reactions 
of two heavy colliding nuclei in experimentally [351 EH] • In accordance with the evolution of the 
colliding system, the whole process of the compound nucleus formation and decay is usually 
divided into three reaction stages, namely the capture process of the colliding system to overcome 
Coulomb barrier, the formation of the compound nucleus to pass over the inner fusion barrier, 
and the de-excitation of the excited compound nucleus against fission [37J. As the first stage 
of synthesizing superheavy nuclei, the accurate calculation of the capture cross sections and the 
analysis of the capture dynamics are very important to estimate the evaporation residue cross 
sections and also affect the competition of the complete fusion and the quasi-fission. Within 
the framework of the ImlQMD model, we analyzed the time evolutions of the relative angular 
momentum and the distance between the centers of the projectile and target in the ^^Ca+^^^U 
reaction at incident energy ii^c.m. =200 MeV as shown in Fig. 10. The relative angular momentum 
is defined as J = {X^^ — Xl^)y/2ilE^ JTi with the reduced mass of the colliding system /i 
and the center of mass distance in the X directions for projectile and target which defines the 
impact parameter at the initial time. The colliding system evolves from the individual nuclei 
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Figure 8: Comparisons of the calculated fusion excitation functions, Hill- Wheeler formula [27] and 
experimental data [251 |22] for the reactions 46'jj_j_46'jj ^^^^^ ^°Ca+^^^Sn. 
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Figure 9: Comparisons of the calculated fusion excitation functions and experimental data for a 
series of reaction systems. 
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Figure 10: The time evolutions of the relative angular momentum and the distance between the 
centers of two colliding nuclei in the reaction ^^Ca+^'^'^U at incident energy 200 MeV. 

to the dinuclear system formation, which is a process of the dissipation of the relative angular 
momentum and the kinetic energy of the relative motion. In Fig. 11 we present a comparison of 
the calculated capture cross sections and the experimental data for a series of reaction systems. 
The experimental data are taken from Refs. [38l [39l HOj HI]. One can see that the calculated 
results are in good agreement with the experimental data. Fig. 12 also shows the calculated 
capture excitation functions for the reactions ^^Ca+^^^Pu and ^^Ca+^^^Cm and compared them 
with available experimental data [12], which have been used to synthesize superheavy elements 
Z=114 and Z=116 in Dubna [13]. Here, we have considered the quadrupole deformation for the 
deformed nucleus at the initial sampling. The capture process of the considered systems takes 
place in the evolution range t=600-800 fm/c. Then, the composite system evolves into either 
separating two fragments with a larger probability (quasi-fission) or maintaining a whole system 
with a smaller probability which leads to the compound nucleus formation. It needs a longer 
evolution time and the larger simulation events to investigate the quasi-fission and the complete 
fusion in the two heavy colliding nuclei. Further works are in progress. 



4 Conclusions 

The shell correction has been further considered in the ImlQMD model. Its influence on the fusion 
excitation functions is analyzed and compared with the experimental data. By using the ImlQMD 
model, the fusion dynamics in heavy-ion collisions near Coulomb barrier is investigated system- 
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Figure 11: Comparisons of the calculated capture excitation functions and experimental data for 
the selected reaction systems. 
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Figure 12: The calculated capture excitation functions for the reactions ^^Ca+^^^Pu and 
^^Ca+^'^^Cm and compared them with available experimental data [12]. 
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atically, such as the fusion barrier, dynamical barrier distribution etc. The calculated excitation 
functions of the fusions and captures are in good agreement with the available experimental data. 
The lowering of the dynamical fusion barrier favors the enhancement of the sub-barrier fusion 
cross sections. The strong shell effect constrains the fusion of two colliding nuclei at near barrier 
energies. 

The physical nature of the heavy-ion fusion reactions near Coulomb barrier is very complicated, 
which is a dynamical process involving the excitation of the colliding system. Semi-classical 
dynamical model can explore its dynamical behavior in a certain extent. There are also some 
interesting works applying the ImlQMD model to investigate the halo nucleus induced reactions, 
the fusion dynamics (capture, quasi-fission and complete fusion) of the synthesis of the superheavy 
nuclei in the massive fusion reactions, and searching other possible ways to produce superheavy 
nuclei etc. For those cases, the structure properties and the longer stability of the initial nucleus 
have to be made. 

5 Acknowledgements 

One of us (Z.Q. Feng) is grateful to Profs. Zhuxia Li, Junqing Li, Werner Scheid, Hans Feldmeier 
and Dr. Ning Wang for fruitful discussions and help, and also thanks the hospitality during 
his stay in GSI. This work was supported by the National Natural Science Foundation of China 
under Grant No. 10475100, the Major state basic research development program under Grant No. 
2007CB815000, and the Helmholtz-DAAD in Germany. 

References 

[1] A.B. Balantekin, N. Takigawa, Rev. Mod. Phys. 70 (1998) 77. 

[2] J. Aichelin, Phys. Rep. 202 (1991) 233. 

[3] E.A. Uehling, G.E. Uhlenbeck, Phys. Rev. 43 (1933) 552. 

[4] L.W. Chen, F.S. Zhang, G.M. Jin, Phys. Rev. C 58 (1998) 2283. 

[5] F.S. Zhang, L.W. Chen, Z.Y. Ming, Z.Y. Zhu, Phys. Rev. C 60 (1999) 064604. 

[6] B.A. Li, CM. Ko, W. Bauer, Int. J. Mod. Phys. E 7 (1998) 147. 



IS 



[7] N. Wang, Z.X. Li, X.Z. Wu, Phys. Rev. C 65 (2002) 064608; N. Wang, X.Z. Wu, Z.X. Li, 
Phys. Rev. C 67 (2003) 024604; N. Wang, Z.X. Li, X.Z. Wu, et al., Phys. Rev. C 69 (2004) 
034608. 



[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 

[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 



T. Neff, H. Feldmeier, K. Langanke, arXiv: |nucl-th/0703030 

Z.Q. Feng, F.S. Zhang, G.M. Jin, X. Huang, Nucl. Phys. A 750 (2005) 232. 

Ch. Hartnack, R.K. Puri, J. Aichehn, et al, Eur. Phys. J. A 1 (1998) 151. 

M. Brack, C. Guet, H.-B. Hakansson, Phys. Rep. 123 (1985) 275. 

E. Chabanat, P. Bonche, P. Haensel, et ah, Nucl. Phys. A 627 (1997) 710; NucL Phys. A 635 
(1998) 231; NucL Phys. A 643 (1998) 441 (erratum). 

Y.X. Zhang, Z.X. Li, Phys. Rev. C 74 (2006) 014602. 

G. Royer, B. Remand, J. Phys. G 10 (1984) 1057. 

Z.Q. Feng, F.S. Zhang, W.F. Li, G.M. Jin, High Ener. Phys. Nucl. Phys., 29 (2005) 41. 

An.R. Leach, Molecular Modelling (Principles and Applications), 1st edition (World Press 
Corporation, Beijing, 1996). 

M. Papa, T. Maruyama, A. Bonasera, Phys. Rev. C 64 (2001) 024612. 

T. Maruyama, A. Ohnishi, H. Horiuchi, Phys. Rev. C 42 (1990) 386. 

N. Wang, Z.X. Li, X.Z. Wu, E.G. Zhao, Mod. Phys. Lett. A 20 (2005) 2619. 

K.A. Brueckner, J.R. Buchler, M.M. Kelly, Phys. Rev. 173 (1968) 944. 

W.D. Myers, W.J. Swiatecki, Phys. Rev. C 62 (2000) 044610. 

K. Siwek-Wilczynska, J. Wilczynski, Phys. Rev. C 64 (2001) 024611. 

Z.Q. Feng, G.M. Jin, F.S. Zhang, et al.. Chin. Phys. Lett. 22 (2005) 3040. 

N. Rowley, G.R. Satchler, P.H. Stelson, Phys. Lett. B 254 (1991) 25. 

C.Y. Wong, Phys. Rev. C 25 (1982) 1460. 

A.M. Stefanini, F. Scarlassara, S. Beghini, et al., Phys. Rev. C 73 (2006) 034606. 

1Q 



[27] D.L. Hill, J.A. Wheeler, Pliys. Rev. 89 (1953) 1102. 

[28] A.M. Stefanini, M. Trotta, L. Corradi, et al, Phys. Rev. C 65 (2002) 034609. 

[29] F. Scarlassara, S. Beghini, G. Montagnoli, et al., Nucl. Phys. A 672 (2000) 99. 

[30] A.M. Borges, C.P. da Silva, D. Pereira, et al., Phys. Rev. C 46 (1992) 2360. 

[31] A.M. Stefanini, L. Corradi, A.M. Vinodkumar, et al., Phys. Rev. C 62 (2000) 014601. 

[32] C.R. Morton, M. Dasgupta, D.J. Hinde, et al., Phys. Rev. Lett. 72 (1994) 4074. 

[33] J.X. Wei, J.R. Leigh, D.J. Hinde, et al., Phys. Rev. Lett. 67 (1991) 3368. 

[34] V.I. Zagrebaev, Phys. Rev. C 67 (2003) R061601. 

[35] S. Hofmann, G. Miinzenberg, Rev. Mod. Phys. 72 (2000) 733; S. Hofmann, Rep. Prog. Phys. 
61 (1998) 639. 

[36] Yu.Ts. Oganessian, J. Phys. G34 (2007) R165; Nucl. Phys. A 787 (2007) 343c. 

[37] Z.Q. Feng, G.M. Jin, F. Fu, J.Q. Li, Nucl. Phys. A 771 (2006) 50; Z.Q. Feng, G.M. Jin, J.Q. 
Li, W. Scheid, Phys. Rev. C 76 (2007) 044606. 

[38] M. Dasgupta, D.J. Hinde, Nucl. Phys. A 734 (2004) 148. 

[39] K. Nishio, H. Ikezoe, Y. Nagame, et al., Phys. Rev. Lett. 93 (2004) 162701. 

[40] E.V. Prokhorova, E.A. Cherepanov, M.G. Itkis, et al., arXiv: |nucl-ex/030902Tj 



[41] W.Q. Shen, J. Albinski, A. Gobbi, et al, Phys. Rev. C 36 (1987) 115. 

[42] M.G. Itkis, J. Ayst5, S. Beghini, et al., Nucl. Phys. A 734 (2004) 136. 

[43] Yu.Ts. Oganessian, V.K. Utyonkov, Yu.V. Lobanov, et al., Phys. Rev. C 62 (2000) R041604; 
Phys. Rev. C 63 (2001) R011301. 



90 



